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We use a perturbative semiclassical trace formula to cal- 
culate the three lowest-order multipole (quadrupole £2, oc- 
tupole £3, and hexadecapole £4) deformations of simple metal 
clusters with 90 < A'' < 550 atoms in their ground states. 
The self-consistent mean field of the valence electrons is mod- 
eled by an axially deformed cavity and the oscillating part of 
the total energy is calculated semiclassically using the short- 
est periodic orbits. The average energy is obtained from a 
liquid-drop model adjusted to the empirical bulk and surface 
properties of the sodium metal. We obtain good qualitative 
agreement with the results of quantum-mechanical calcula- 
tions using Strutinsky's shell-correction method. 

PACS numbers: 03.65. Sq, 05.30.Fk, 31.15.Ew, 71.10.Ca 

Free clusters made of simple metal atoms exhibit a pro- 
nounced electronic shell structure Although the 
detailed experimental information obtained, e.g., from 
photo-excitation measurements can only be understood 
if the ionic structure is taken into account ^ , the quali- 
tative features of the electronic shell structure can be well 
described, for not too small systems, by phenomcnologi- 
cal deformed shell-model potentials . Self-consistent 
density functional calculations in the framework of a de- 
formed jellium model |^ have revealed that the cluster 
ground-state shapes can be well characterized in terms of 
the lowest three multipole orders £2 (quadrupole) , £3 (oc- 
tupole) , and £4 (hexadecapole) . Since such self-consistent 
calculations are quite time consuming computationally, it 
is often more efficient to resort to simpler methods, such 
as the shell-correction method introduced by Strutinsky 
in nuclear physics in particular, if more shape degrees 
of freedom are to be investigated 0. 

An even more economical approach is the semiclassi- 
cal periodic orbit theory (POT) (see, e.g., Ref. jl^ for 
a general introduction), in which quantum oscillations in 
the level density or other observables can be described in 
terms of the leading shortest periodic orbits of the corre- 
sponding classical system through so-called trace formu- 
lae 12|. This method has been used for quadrupole- 
deformed clusters in a Nilsson-type model ||l3|| and, more 
recently, using cavities with axial £2, 63, and £4 deforma- 
tions The approximation of the self-consistent 
mean field of the valence electrons by a cavity with re- 
flecting walls has received strong support from the quan- 
titative explanation iQ of the experimental magic num- 
bers found in connection with the electronic supershells 
pa in terms of the trace formula of the spherical cavity 



p2[ . The validity of the cavity model has also been con- 
firmed by calculations with more realistic Woods-Saxon 
type potentials and by selfconsistent Kohn-Sham cal- 
culations pql in the spherical jellium model. 

In Re f. ||15| , a perturbative trace formula derived by 
Creagh [|19[| has been used for axially deformed cavities 
with small multipole deformations £2, £3, or £4 and found 
to reproduce the quantum-mechanical results very well 
for moderate values of the deformation parameters. The 
main virtue of this approach is that the perturbed orbits 
need not be known. The trace formula is given in terms of 
the periodic orbits of the unperturbed system (i.e., here, 
the spherical cavity), their unperturbed amplitudes be- 
ing multiplied by a modulation factor that contains the 
perturbations of their action integrals. The periodic or- 
bits of the spherical cavity were discussed first by Balian 
and Bloch |l2[ and characterized by their winding num- 
ber t (number of revolutions around the centre) and the 
number p of reflections from the boundary. 

In the present work we extend the approach developed 
in Ref. The deformed mean field of the electrons 

is parametrized by an axial cavity whose surface is de- 
scribed in terms of the Legendre polynomials Pj by 
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R{e) = R(l + eo + Y^ £jFj (cos( 



(1) 



where R — r^N^^^ is the radius of the spherical cavity 
and £0 is used to conserve its volume. In contrast to 
Ref. , we include here simultaneously all three defor- 
mation parameters (def) — (£2, £3, £4) and minimize for 
each cluster size N the total energy 

Etot{N, def) = Eldm{N, def) + SE{N, def) , (2) 

in order to determine its ground-state deformation. Here 
Eldm is the average total energy obtained in the liq- 
uid drop model (LDM) and 5E is the shell-correction 
energy. For the LDM model parameters and the radius 
constant we use the same values as in Ref. Q . We eval- 
uate 5E both quantum-mechanically using the Strutin- 
sky method, and semiclassically by the perturbative trace 
formula discussed in Ref. . 

The quantum-mechanical energy spectrum {Ei{def)} 
of the spheroidal cavities with the shapes defined by (|l|) 
was obtained by a method developed in Ref. In 
terms of the Ei{def), the exact (quantum- mechanical) 
shell-correction energy is then defined as usual M by 
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SEqm ^Y.^^- ^ 9{E) dE , (3) 
i=i ■'0 

whereby 'g{E) and A are given further below. 

We refer to Rcf. for the definition of the mod- 
ulation factor included in the perturbative semiclassi- 
cal trace formula for the oscillating part Sgsc{E) of the 
level density. When admitting several types of defor- 
mations ej simultaneously, the corresponding first-order 
changes of the classical actions ASj have to be added 
up in the exponent of the integrand for the modulation 
factor. The semiclassical shell-correction energy SEgc is 
then obtained using the relation 

6Esc= / E[g{E)+6gsc{E)]dE- / Eg{E)dE, (4) 
Jo Jo 

where the Fermi energies A and A are determined itera- 
tively by conserving the particle number: 

N = [ [g{E)+6gsc{E)]dE^ [ g{E) dE . (5) 
Jo Jo 

Here g{E) is the average level density found analytically 
from its Weyl expansion (see, e.g., Ref. [0) or by a nu- 
merical Strutinsky averaging |^ of the spectrum {Ei}. 
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FIG. 1. Deformation energy versus atom number N. 
Circles: quantum-mechanical with spherical shapes; crosses: 
quantum-mechanical, deformed. Solid line: semiclassical. For 
the latter two, the energy was minimized for each A'^ with re- 
spect to the three multipole deformations e^, es, and 64. 

In Fig. 1^ we present deformation energies E^ciidef) — 
Etot{def) — EL]jM{0) versus the particle number N in the 
region of 90 < < 200. The circles show the quantum- 
mechanical Strutinsky results for spherical shapes and 
exhibit three typical shell closures at the magic numbers 
N = 92, 138, and 186. The crosses show the quantum 
results obtained with deformed shapes, minimizing Etot 
for each N with respect to the three multipole deforma- 
tions £2, £3, and £4. The solid line is the result of the 



semiclassical calculation including periodic orbits with 
up to tm = 3 repetitions and Pm = 30 reflections. We 
see how the inclusion of deformations reduces the total 
energy drastically, except for the narrow regions around 
the magic numbers. We obtain a very good agreement of 
our semiclassical results with the quantum results, show- 
ing that the perturbative trace formula works reasonably 
well even for the largest deformations. 
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FIG. 2. Ground-state quadrupole (£2), octupole (£3) and 
hexadecapole (£4) deformations of Na clusters versus number 
of atoms A*'. Crosses: quantum-mechanical results. Lines give 
the semiclassical results with different numbers of included 
periodic orbits. Solid lines: tm ~ 3, pm = 30; dashed lines: 
tm = 1, Pm = 10 (see text). 

In Fig. H, we present the three ground-state deforma- 
tions £2, £3, and £4, obtained by the energy minimiza- 
tion procedure, as functions of N . The crosses are the 
quantum-mechanical results, and the solid and dashed 
lines are the semiclassical results. The agreement be- 
tween semiclassics and quantum mechanics is very sat- 
isfactory in view of the fact that we cannot expect the 
semiclassical approximation to be exact. In particular, 
even the largest deformations agree within about 10% - 
20%, in spite of the fact that for these deformations the 
argument of the modulation factor in the trace formula 
is larger than unity. A particularly nice result is the 
correct reproduction of the systematics of octupole de- 
formations £3 which only occur in the lowest part of each 
shell, as already noticed in Ref. |^]. In addition, the sys- 
tematic sign change in the hexadecapole deformations £4 
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from positive to negative within each shell is reproduced 
qualitatively, although the negative values turn out to 
be systematically too small in the semiclassical results. 
This defect may be connected to the fact that for larger 
negative hexadecapole deformations the cluster shape be- 
comes a multi-valued function in the cylindrical coordi- 
nates. Some minor changes of the deformation system- 
atics are expected when non-axial shapes are used 
although they can hardly be confirmed experimentally. 

The convergence of the semiclassical results with re- 
spect to the lengths of the included periodic orbits is 
demonstrated by the dashed line in Fig. ^ for which 
(trmPm) = (IjlO) was used, whereas the solid line was 
obtained for {tm,Pm) = (3,30) as in Fig. ^ We clearly 
see that a slight improvement is obtained by including the 
second {t = 2) and third {t = 3) repetitions (harmonics) 
and the corresponding numbers of reflections. 

With the same choice of maximum orbit lengths, wc 
have calculated the deformation energies and ground- 
state deformations of clusters with up to = 550 atoms. 
The results are presented in Figs. and ^. Here the same 
systematics are found as in Fig. above, although they 
become less distinct for the larger clusters where the su- 
pershells start to decrease in amplitude. 
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FIG. 3. Deformation energy (in eV) versus atom number 
iV. Crosses: quantum-mechanical with spherical shapes; solid 
line: optimized semiclassical results {tm = 3, pm ~ 30). 

Owing to the efficiency of our method, the numerical 
calculation of multidimensional deformation energies be- 
comes quite economical. As examples, we show in Figs. 
^ and 1^ two-dimensional deformation energy surfaces of 
Nai62 and Naigo- The deformation energy Edef is plot- 
ted versus quadrupole (62) and hexadecapole (£4) defor- 
mations. The upper parts give the quantum-mechanical 
results and the lower parts the semiclassical ones. We see 
that our semiclassical method captures the correct over- 
all behaviour of the topology of the quantum-mechanical 
deformation energy landscape, thus being able to repro- 
duce the most prominent isomeric minima correctly. 



We have purposely limited our investigations to clus- 
ter sizes A <; 90. Several reasons limit our approach to 
such particle numbers. First, the selfconsistent total po- 
tential (which includes the Coulomb repulsion between 
the electrons) does not have very steep walls for particle 
numbers A < 60 - 80 (cf. |l|,|l|,|l^), so that the cavity 
model cannot be justified. Second, the ionic structure ef- 
fects become more important for smaller clusters, so that 
any smooth shell-model potential becomes a bad approx- 
imation. Third, the absolute values of the ground-state 
deformations tend to increase with decreasing particle 
number , so that the perturbative approach will break 
down for small clusters. 
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FIG. 4. Ground-state deformations for larger clusters with 
200 < N < 550, obtained semiclassically (tm = 3, pm = 30). 



In summary, we have shown that the perturbative 
semiclassical POT is an efficient and reliable tool for the 
calculation of ground-state deformations of metal clus- 
ters, or of any other system of fermions whose mean field 
can be approximated by a cavity with refiecting walls. 
This holds also for atomic nuclei ||2^ , apart from modifi- 
cations due to the spin-orbit interaction whose inclusion 
in the POT is still an object of actual research p3| , p4| . 
Our method can easily be extended to describe charged 
clusters by using the relevant liquid drop model devel- 
oped, e.g., in Ref. ||25|] . 

We are grateful to S. Fraucndorf for his interest and 
encouraging discussions at early stages of this work. P.M. 
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FIG. 5. Deformation energy Edef of Nai62 versus deforma- 
tions £2 and £4. Upper part: quantum-mechanical, lower part: 
semiclassical result {tm = 3, Pm ~ 30). 
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